Long Range Correlations and the Soft Ridge in Relativistic Nuclear Collisions 
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Relativistic Heavy Ion Collider experiments exhibit correlations peaked in relative azimuthal angle 
and extended in rapidity. Called the ridge, this peak occurs both with and without a jet trigger. 
We argue that the untriggered ridge arises when particles formed by flux tubes in an early Glasma 
stage later manifest transverse flow. Combining a blast wave model of flow fixed by single-particle 
spectra with a simple description of the Glasma, we find excellent agreement with current data. 
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Correlation measurements of high transverse momen- 
tum particles provided the first striking experimental ev- 
idence of jet quenching. Further studies of the correla- 
tion of high p t particles with more typical particles re- 
veal a complex pattern of correlations as a function of 
relative pseudorapidity rj = r\\ — r\i and azimuthal angle 
4> = <pi — 4>2- In particular, a "hard ridge" of enhanced 
correlation is observed near r\ = (f> = that is narrow in 
4> and broad in rj pQ. Interestingly, the STAR collabo- 
ration reports a similar ridge in correlations of particles 
of any p t , i.e., without a jet trigger [21 El H] . As with 
the hard ridge, the width of the untriggered soft ridge is 
broad in rj and narrow in 4>. Preliminary PHOBOS data 
suggests that the hard ridge and possibly the soft ridge 
may extend over the broad range —4 < r/ < 2 [S]. 

We argue that the soft ridge is a consequence of early- 
stage rapidity correlations in concert with late-stage 
transverse flow. Correlations over several rapidity units 
can only originate at the earliest stages of an ion col- 
lision when the first partons are produced [B]. Hydro- 
dynamics and other later-stage effects can modify these 
correlations, but are limited to a horizon of ~ 1 — 2 rapid- 
ity units. Analogous to super-horizon fluctuations in the 
cosmos, these long range correlations can therefore reveal 
the 'little bang' in each nuclear collision at its birth. 

Almost instantaneously after a collision of two nuclei, 
Color Glass Condensate (CGC) theory predicts that the 
transverse fields of each nucleus are transformed into lon- 
gitudinal fields that are approximately uniform in rapid- 
ity [3 [8]. The fields, which are random over transverse 
distances larger than Qj 1 where Q s ~ 1 — 2 GeV is the 
saturation scale, comprise a series of flux tubes. Long 
range rapidity correlations arise because particles from 
the same flux tube start at nearly the same transverse 
position, regardless of rapidity. 

Pressure builds as the flux tubes fragment to form par- 
ticles. The resulting transverse flow modifies these long 
range correlations by focusing particles into a narrow re- 
gion in (j>. Suppose that the transverse fluid velocity has 
the Hubble- like form 7 t v t = Ar t . A fluid cell a distance 



r t from the center of the collision volume will then have 
a mean speed v t - When this cell freezes out, it releases 
particles into an opening angle <fi ~ Vth/vt oc (Ar t ) _1 , 
for a thermal velocity v t h ~ 1. This modifies the cor- 
relations, since particles near the center of the collision 
volume have a large opening angle, while those from a 
larger r t have a smaller <f>. Voloshin has long stressed the 
connection between flow and pt correlations [3] . 

In this letter we focus on the centrality dependence 
of the height and azimuthal width of the near side peak 
of the soft ridge. We will see that this dependence can 
be explained using CGC-Glasma scaling arguments com- 
bined with blast wave calculations. However, we remark 
that the measured correlation function is not a plateau 
in rapidity, but a broad structure of width ^1 — 2 units, 
perhaps with tails extending higher in rapidity |H |5J. 
We will address this rapidity dependence elsewhere, since 
that analysis requires quantum corrections to the Glasma 
[6j as well as viscous corrections to the hydrodynamic 
treatment [TU] . We will also not discuss the jet triggered 
data, which would require a description of the passage of 
the jet through the high-density environment produced 
by the nuclear collision [TTJ [T2] . 

Flux tubes arise naturally in descriptions of high en- 
ergy collisions [T3] . The Glasma description incorporates 
many of these features, and in a high density environment 
such as is produced by the collisions of nuclei, allows for 
a systematic weak coupling computation. The contribu- 
tion of flux tubes to long range correlations is studied 
in the Glasma formulation in ref. [5J. We imagine the 
Glasma to be filled with flux tubes of large longitudinal 
extent but small transverse size ~ Qj • Each flux tube 
yields a multiplicity of ~ a s {Q s )^ 1 gluons. The number 
of flux tubes is proportional to the transverse area R\ 
divided by the area per flux tube, Q~ 2 . The rapidity 
density of gluons is therefore 

dN/dy~a s - l Q 2 s R 2 A . (1) 

The number of final hadrons scales similarly [14 . 

We characterize correlations in the Glasma and later 
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in the evolution using the spatial correlation function 

c(xi,x 2 ) = ti 2 (xi,x 2 ) - ni(xi)ni(x 2 ), (2) 

where n\ and n 2 are the single and pair densities. In the 
absence of correlations, n 2 (xi,x 2 ) — > rii(xi)ni(x 2 ) so 
that c vanishes. The integral n 2 over both positions gives 
the number of pairs averaged over events (N(N — 1)). 
When correlations are negligible, the integral of c van- 
ishes - as it must - because N follows Poisson statistics 
and, therefore, (N(N-l)) -> (N) 2 . 

We take pairs from the same flux tube as correlated 
and neglect correlations between tubes. Furthermore, 
we assume that correlations are independent of rapid- 
ity. The correlation function then depends only on the 
relative transverse position r t = t — t as well as the 
average R t = (ri. t + r i, t)/2- The correlation length in r t 
is roughly the flux tube size ~ Q^ 1 , while the correlation 
length in R t of order of the transverse system size Ra- 
For Qj 1 <§; Ra we take the correlation function to be 
point- like in r t and broad in R t , writing 



c(xi,x 2 ) =U5(r t )p FT (R t ). 



Here, p FT (Rj) describes the transverse distribution of 
flux tubes in the collision volume, which we assume fol- 
lows the thickness function of the colliding nuclei 



a (R)- 2{N)2 fl R * 

pft( r *) - -^jt y r\ 



(4) 



for R t < Ra, and zero otherwise. Integrating both sides 
of eq. ^ with respect to r t and R t , we find 



(N) 2 K 



cd 3 Xl d 3 x 2 = (N 2 ) - (N) 2 - (N). (5) 



To see how 1Z depends on Q s , think of each flux tube as a 
source that produces particles with a mean multiplicity fi 
and variance a 2 . For K flux tubes, the mean multiplicity 
is ixK and the variance is a 2 K. If K fluctuates from 
event to event, then the mean multiplicity is p(K) and 
the variance is a 2 (K) + ^ 2 {{K 2 ) - (K) 2 ). Therefore 



K = 



1 



a 2 (K) 



(K 2 ) (KY< 
(K) 2 



(6) 



Particle production from a flux tube is a Poisson process, 
since the flux tube is a coherent state. It follows that 
(j 2 = fj,, so that the first contribution vanishes. For large 
K, the second term is oc (K)^ 1 . 

We combine these results to obtain a scaling relation 
for the integrated strength of correlations in the Glasma 



Hoc (K)- 1 = (Q S R)~ 



(7) 



a result supported by momentum-space calculations in 
[6J. In contrast, the mean multiplicity in a rapidity inter- 
val scales as a s (Q s )~ 1 Q 2 R 2 ; see eq. ([TJ. This difference 



will prove significant later. We comment that ^ and 
similar CGC relations may not quantitatively describe 
pp or peripheral collisions at the energies studied here, 
although phcnomenological string models may apply. 

We now turn to discuss the impact of these long range 
correlations on the final-state particle correlations. As 
the partons emitted from these flux tubes locally equi- 
librate, transverse flow builds. To describe the effect of 
thermalization and flow on the pair correlation function 
at freeze out, we generalize the common blast- wave model 
US [igm QBE]]. To begin, recall that the Cooper-Frye 
single-particle distribution is 



Pl (p) = dN/dyd 2 Pt = J /(x,p)dT, 



(8) 



where /(x, p) = (2ir)~ 3 exp{— p^u^/T} is the Boltzmann 
phase-space density for a temperature T and fluid four- 
velocity Up, and dT = p^da^ is the element of flux 
through the four dimensional freeze out surface along 
which particle interactions effectively cease. We assume 
that freeze out occurs at a proper time tf, so that 
(3) P^da^ — Tprrit cosh(y — rj)drjd 2 r t , where rj = (l/2)ln((t + 
z)/(t — z)) is the spatial rapidity. We follow ref. [TS] and 
write the four velocity of the longitudinal-boost invari- 



ant blast wave as u M = 74 (cosh 77, v t ,sinh?7), where Vt 
is the transverse velocity and j t = (1 — v 2 )^ 1 ^ 2 . The 
phase space density is then / oc exp{— 7tWt cosh(y — 
T])/T} exp{7 t v t • pt/T}. We take the transverse velocity 
to be 7tV t w Ar t , a widely-used ansatz that adequately 
describes much of SPS and RHIC data. The calculation 
of ^ is standard and follows [15] . 

To exhibit the effect of flow on particle correlations, 
we use the momentum-space correlation function 



KPl>Pa) = P2(Pl,P2) - Pl(Pl)Pl(P2) 



(9) 



where p 2 (pi,p 2 ) = dN j dy\d 2 p t \dy2d 2 pti is the pair dis- 
tribution. Generalizing (JsJ) , we write 

I \ f I n/( x 1;Pi) /( x 2,P2) jt, / 1n x 

KPi>P2j = / c(x 1 ,x 2 ) — r — — dT 1 dT 2 . (10) 

J ni(xi) ni(x 2 ) 

We identify c(xi,x 2 ) at freeze out with ([3]), a form that 
describes the system at its formation. This identifica- 
tion omits the effects of diffusion described in ref. |10j . 
This omission is reasonable only as long as we restrict 
our attention to the long range correlations with pairs 
separated by \r)i — 372 1 > 1. 

STAR measures the characteristics of the untagged 
near-side ridge as functions of the centrality at 62 and 
200 GeV for Au+Au @]. While they focus on the re- 
gion — 1 < 77 < 1 where short and long range correla- 
tion phenomena arc both present, it is instructive to see 
which aspects of the data can be explained by a purely 
long range model. To facilitate or comparison, we vi- 
sualize the STAR analysis as consisting of the follow- 
ing steps. First, a correlated two particle distribution 
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of "sibling" particles p S i is measured. This quantity 
is essentially our p^ integrated over the magnitudes of 
each particle's pt as well as the average azimuthal angle 
$ = {(j>i + 4>%)l2 and pseudorapidity r\ a = (r\\ + rj 2 )/2. 
The resulting density depends only on the relative quan- 
tities <p — 0i — <f>2 and 77 = 771 — 772- Second, an uncorre- 
cted pair distribution p re f is obtained from mixed events 
and Ap/yfp = (p s a, — pref)/ \j Pref is constructed. Next, 
a rapidity-independent function a + b cos + c cos 20 is 
subtracted to remove backgrounds as well as elliptical 
flow and momentum conservation effects. Finally, the 
corrected (77, 0) distribution is subjected to a multicom- 
ponent fit to extract the attributes of the near side peak. 
In practice, these steps are performed simultaneously. 

To confront the STAR measurements, we calculate 
Ap(j], 0) by integrating ( 10 I over all momenta. Similarly, 



we compute p re f(r],<fi) by integrating Pi(pi)pi(p2)- The 
integrations are straightforward if we take the 7 factors 
as constants evaluated at R t — Ra- We obtain 



Ap/^p^j = KdN/dy F(0), 
where f Q * F(0)<i0 = 1. The factor 

TZdN/dy = HiasiQs)- 1 



(11) 



(12) 



follows from the Glasma relations ([T]) and Q, where k is 
an energy independent constant to be determined from 
data. The angular distribution F(<f>) depends only on 
blast wave parameters "fm/T and v s . 

We specify the centrality dependence of the correlation 
function using the velocity and temperature fit from sin- 
gle particle spectra at 200 GeV in ref . [16] . The computed 
Ap/^fp is shown as the dashed line in the top panel of 
fig. 1. We fix k in (12 1 and (111 to agree with the mag- 
nitude of the 200 GeV data. We define the height of 
the near side peak as the difference between Ap/^fp at 
0=0 and 7r. The dashed line in the top panel is the 
blast wave result F(<p) without the CGC scaling. It fol- 
lows the basic trend of the data rather well, given that 
the only parameters that vary with centrality - v and T 
— are fit elsewhere |16| . The uncertainty in v and T im- 
plies the shaded bands in fig. 1. To compute Ap/^/p at 
62 GeV, we follow [TB] and reduce the velocities by 5% 
and the temperatures by 10% by uniform scale factors. 
The dashed curve in the bottom panel of fig. 1 is well 
above the data, but agrees roughly in shape. 

The very existence of long range correlations implies 
strong correlations in the initial state, as predicted by 
Color Glass Condensate theory [BJ. That said, it is in- 
teresting to see how (12 1 influences the systematics of 
the soft ridge. This modifies the centrality dependence 
of Ap/y/p by introducing a logarithmic dependence on 
Q s , which in turn depends on N part . It is important 
to note that although the dashed curve represents the 
blast wave without the a s scaling, the correlation func- 
tion in the blast wave integrals enforces restraints from 
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FIG. 1: Height of the near side peak vs. centrality for 200 
GeV (top) and 62 GeV (bottom). Preliminary STAR data 
is from [I]. Bands indicate the uncertainty of the blast wave 
parameters T and v. 



the Glasma picture. The upper solid curve in fig. 1 com- 
bines this a s (Npart) -1 dependence from ref. [2] with the 
blast wave behavior. Agreement with data is impressive 
given that Q s (Np art ) is obtained fit in |14j . 

We deduce Ap/^/p at 62 GeV using the Q s dependence 
discussed earlier together with the relevant blast wave v 



and T. With the proportionality constant in ( 12 1 fixed by 
the 200 GeV data, there are no further free parameters 
to adjust. The solid curve in the lower panel on fig. 1 
is in good accord with the data for the expected drop of 
Q 2 S by ~ 1/2 relative to the 200 GeV value [14]. Observe 
that most of the change from 62 to 200 GeV owes to the 
ctjHQs) dependence, since the change in the blast wave 
parameters is small [15] , 

We compare our calculations to the measured az- 
imuthal width of the near side peak in fig. 2. To simulate 
the experimental fit procedure, we obtain this width by 
fitting a gaussian plus a constant offset to the computed 
distribution in the near-side interval — 7r/2 < < 7r/2. 
The uncertainty band in fig. 2 indicates the impact of 
changing the near-side interval by ±20%. Once again, 
the agreement is surprisingly good given the simplicity of 
the model. The calculated angular width does not change 
in this energy range, since the normalization does not af- 
fect F((f>). Note that ref. [B] includes flow by boosting 
a Glasma source. Their computed width exceeds data 
because they omit the pt dependence of the boost and do 
not simulate the experimental fit procedure. 

In summary, correlation measurements show a ridge 
that is narrow in and broad in 77, perhaps extend- 
ing several units in rapidity. Such long-range rapidity 
correlations can only be caused by super-horizon fluctu- 
ations at sub-fermi time scales. While flux tubes may 
also enhance forward-backward correlations [20]. These 
measurements therefore provide an image of the particle 
production process at the sub-fermi scale, which can be 
corroborated e.g., by forward-backward correlation mea- 
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FIG. 2: Azimuthal width of the near side peak vs. centrality. 
Preliminary STAR data is from [4]. The curve is obtained 
for by fitting a gaussian to the computed distribution in the 
range —n/2 < <f> < n/2. The band shows the sensitivity of 
the result to a 20% change in this range. 

surements |20j . For simplicity we have focused on the 
soft ridge, but similar considerations may also apply to 
jet-tagged measurements [T^]. The height and azimuthal 
width measured near midrapidity in ref. [1] are consistent 
with flux tubes of large longitudinal extent formed early 
in nuclear collisions. Correlations predicted by Color 
Glass Condensate theory combined with transverse flow 
provide a remarkably good description of the near side 
ridge. In particular, agreement with 62 GeV data fol- 
lows mainly from the CGC-Glasma prediction ([7]), the 
only free parameter being an overall constant fixed at 
200 GeV. This agreement is surprising, since our model 
only includes long range correlations. The rapidity de- 
pendence for \r)\ < 1 — 2 requires a more detailed hydro- 
dynamic description [TP] . 
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